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Abstract. 

We present a survey of interplanetary (IP) shocks using WIND and ACE satellite data 
from January 1995 to December 2013 to study how IP shock geoeffectiveness is controlled 
by IP shock impact angles. A shock list covering one and a half solar cycle is compiled. 

The yearly number of IP shocks is found to correlate well with the monthly sunspot num¬ 
ber. We use data from SuperMAG, a large chain with more than 300 geomagnetic sta¬ 
tions, to study geoeffectiveness triggered by IP shocks. The SuperMAG SML index, an 
enhanced version of the familiar AL index, is used in our statistical analysis. The jumps 
of the SML index triggered by IP shock impacts on the Earth’s magnetosphere is inves¬ 
tigated in terms of IP shock orientation and speed. We find that, in general, strong (high 
speed) and almost frontal (small impact angle) shocks are more geoeffective than inclined 
shocks with low speed. The strongest correlation (correlation coefficient R = 0.70) oc¬ 
curs for fixed IP shock speed and varying the IP shock impact angle. We attribute this 
result, predicted previously with simulations, to the fact that frontal shocks compress 
the magnetosphere symmetrically from all sides, which is a favorable condition for the 
release of magnetic energy stored in the magnetotail, which in turn can produce mod¬ 
erate to strong auroral substorms, which are then observed by ground based magnetome¬ 
ters. 


1. Introduction 

Interplanetary (IP) shocks occur throughout the helio¬ 
sphere as a result of the interaction of solar disturbances 
with the solar wind [Burlaga, 1971; Richter et al., 1985]. 
As IP shocks interact with the Earth’s magnetosphere, they 
cause disturbances that can be seen throughout the mag¬ 
netosphere. Some of these disturbances can have implica¬ 
tions in several sectors of both magnetosphere and iono¬ 
sphere [Oliveira, 2014], for example, SSCs (storm sudden 
commencements), geomagnetic storms, auroral substorms, 
and GICs (ground-induced currents). GICs may impact 
power grids, causing electric power disruptions due to equip¬ 
ment damage [Bolduc, 2002; Kappenman, 2010], and inter¬ 
ruption of commercial activities leading to severe economic 
losses [Schrijver et al., 2014]. 

At 1 AU, most IP shocks are fast shocks. Although slow 
shocks occur close to the Sun, a few slow shocks may be 
observed at Earth’s orbit [ Chao and Olhert, 1970; Whang 
et al, 1996]. IP shocks are then classified as forward and 
reverse. Forward shocks propagate away from the Sun, and 
reverse shocks propagate towards the Sun in the solar wind 
frame. Since the solar wind speed is almost always super- 
magnetosonic, all shocks propagate away from the Sun in the 
Earth’s frame. IP shocks may be further classified by their 
strength in terms of Mach numbers and the compression 
ratio, the ratio of downstream to upstream plasma densi¬ 
ties. Among other parameters, the shock normal is another 
important feature of IP shocks, because shock normal orien¬ 
tations determine how IP shocks propagate throughout the 
heliosphere. Normals of most IP shocks generated by ICMEs 
(the interplanetary manifestation of coronal mass ejections 
at the Sun) at 1 AU are concentrated near the Sun-Earth 
line [Richter et al., 1985]. However, shocks driven by coro¬ 
tating interaction regions (CIRs), as a result of the slow solar 
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wind compression by a fast stream, have normals inclined in 
relation to the Sun-Earth line [see, e.g. Siscoe, 1976; Pizzo, 
1991, and references therein]. For CIR-driven shocks, the 
normal angles in the azimuthal direction in relation to the 
solar coordinate system are generally equal or larger than 
the inclination angle [Siscoe, 1976; Pizzo, 1991]. Calcula¬ 
tions of IP shock normals are very sensitive to upstream 
and downstream plasma parameters, and using data from 
more than one spacecraft is thought to improve shock nor¬ 
mal determinations [Russell et al., 2000]. 

Here our primary concern is to study the influence of 
IP shock normal orientations on IP shock geoeffectiveness. 
Similar studies have been done in the past, but with differ¬ 
ent space parameters. For example, Jurac et al. [2002] found 
that the angle between the shock normal and the upstream 
magnetic held plays an important role in the shock geoef¬ 
fectiveness. Their main result was the finding that quasi¬ 
perpendicular shocks were more geoeffective than quasi¬ 
parallel shocks. Takeuchi et al. [2002] showed that an IP 
shock highly inclined in relation to the equatorial plane led 
to an unusual SSC rise time (~ 30 min). They attributed 
this result to the fact that the IP shock observed took a 
longer time to sweep over what they called the “geoeffec- 
tive magnetopause”. By performing global MHD simula¬ 
tions, Guo et al. [2005] showed that two similar IP shocks 
with different shock normal orientations, namely a frontal 
shock and a highly inclined shock whose shock normals lay 
in the equatorial plane, interacted with the Earth’s mag¬ 
netosphere leading to different conclusions. They reported 
that the inclined case took longer than the frontal case to 
reach fairly similar final quasi-steady states. Later on, Wang 
et al. [2006], using ACE and WIND satellite data from 1995 
to 2004, reported that, in a survey of nearly 300 fast forward 
IP shocks, 75% of them were followed by SSCs observed on 
the ground. They also found that the shock impact angle 
plays an important role in determining the SSC rise time, 
as previously suggested by Takeuchi et al. [2002]. When the 
shock speed (shock strength) was fixed, the more parallel 
the shock normal with the Sun-Earth line, the smaller the 
rise time. The same occurred when they fixed the shock 
inclination and changed the shock speed. The faster the 
shock, the shorter the SSC rise time. 
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More recently, the geoeffectiveness of IP shocks controlled 
by IP shock impact angles was studied by Oliveira and 
Raeder [2014] through global MHD simulations. Using the 
OpenGGCM (Open Global Geospace Circulation Model) 
MHD code [Raeder, 2003], they showed that similar IP 
shocks with different IP shock impact angles may lead to 
different IP shock geoeffectiveness. They simulated three 
different IP shocks, where two had shock normals inclined 
in relation to the Sun-Earth line in the meridian plane, while 
the second had a Mach number as twice as the first. Both 
shocks were oblique, i.e., their shock normals were at angles 
close to 45° with the upstream magnetic field in the shock 
frame of reference. Finally, the impact of a third shock, 
namely a frontal and perpendicular shock with the same 
Mach number as the first shock, with the Earth’s magne¬ 
tosphere was simulated. Oliveira and Raeder [2014] found 
that the third shock was much more geoeffective than the 
other two because the shock was frontal, and the magne¬ 
tosphere was compressed symmetrically on both north and 
south sides. This compression led then to the triggering of 
a strong auroral substorm not seen in the other cases. 

The goal of this work is to confirm the role of IP shock im¬ 
pact angles in the IP shock geoeffectiveness using satellite 
and geomagnetic activity data. The geomagnetic activity 
triggered by such interactions is then analyzed using an en¬ 
hanced version of the geomagnetic AL index. The data and 
methodology used are discussed in section 2. We present our 
statistical results in section 3, and in section 4 we summarize 
and discuss our results. 

2. Data and methodology 

2.1. Data 

We screen solar wind plasma and field data at 1 AU to find 
fast forward IP shock events. In order to do so, we use two 
different spacecraft close to the equatorial plane: WIND, 
with data from 1995 up to 2013, and ACE (Advanced Com¬ 
position Explorer), with data from 1998 also up to 2013. The 
WIND data were obtained from the Solar Wind Experiment 
(SWE) instrument [Ogilvie et al., 1995], and the Magnetic 
Field Investigation (MFI) instrument [Lepping et al, 1995], 
both in 93-second time resolution. The ACE data were ob¬ 
tained from the Solar Wind Electron Proton Alpha Moni¬ 
tor (SWEPAM) instrument [ McComas et al., 1998] and the 
ACE Magnetic Field Experiment (MAG) instrument [ Smith 
et al., 1998], both with 64-second time resolution. All data 
were downloaded from the CDAWeb interface located at 
http://cdaweb.gsfc.nasa.gov. All these data were used to 
compile a list of 461 fast forward IP shock events that can 
be found in the supplementary material. 

The monthly averaged sunspot number (SSN) data 
were compiled by the Solar Influence Data Analy¬ 
sis Center (SIDC). This list was downloaded from 
http://side.oma.be/silso/datafiles. 

It is well established that substorm activity may be trig¬ 
gered by IP shocks [Kokubun et al, 1977; Akasofu and Chao, 
1980; Zhou and Tsurutani, 2001; Yue et al., 2010; Tsuru- 
tani et al, 2014; Oliveira and Raeder, 2014], and that AL 
appears to be the preferred index to quantify the strength 
of auroral activity [Mayaud, 1980]. The AE index, the auro¬ 
ral electroject index, was first suggested by Davis and Sug- 
iura [1966] and has been heavily used by magnetospheric 
physicists since then. However, as pointed out by Davis 
and Sugiura [1966] themselves and reviewed by Rostoker 
[1972], the indices AU, AL, and AE = AU - AL, available 
at the World Data Center (WDC) in Kyoto, Japan, web¬ 
site (http://wdc.kugi.kyoto-u.ac.jp/aeasy/index.html.), are 
limited because of the relatively low number of ground sta¬ 
tions used to define these indices which happens to be 12. 


Therefore it is clear that sometimes strong auroral events 
are underestimated because there are no ground stations 
under the auroral bulge contributing to the construction of 
these indices during some strong auroral substorm events 
[Gjerloev et al., 2004]. As an alternative to alleviate this 
deficiency, SuperMAG, a large worldwide collaboration in¬ 
volving more than 300 ground-based magnetometers, was 
formed [Gjerloev, 2009]. Because the AU, AL, and AE in¬ 
dices are recognized as official indices by IAGA, SuperMAG 
defines SMU as the SuperMAG measurement of the maxi¬ 
mum eastward auroral electrojet strength (upper envelope 
of N-component measured by stations between 40° and 80° 
magnetic north), SML as the SuperMAG measurement of 
the minimum westward auroral electrojet strength (lower 
envelope of N-component measured by stations between 40° 
and 80° magnetic north), and SME = SMU - SML as the 
SuperMAG measurement of the auroral electroject index de¬ 
fined as the distance between the last two indices [Newell 
and Gjerloev, 2011]. 

An example of an auroral substorm event observed by 
different numbers of IAGA and SuperMAG stations is rep¬ 
resented by Figure 1 in Newell and Gjerloev [2011]. In their 
event, it is shown by Polar UVI imagery that the expansion 
of the auroral bulge does not pass over any AE ground sta¬ 
tions, but rather passes over more than half dozen of the 
SME ground stations. This auroral substorm was underes¬ 
timated by the AE stations, as shown by Newell and Gjer¬ 
loev [2011] in their Figure 2. The AL stations did not de¬ 
tect a substorm event; however, the SML stations recorded 
a substorm onset 37 seconds after the onset registered by 
Polar UVI observations. Therefore it is important to men¬ 
tion that AE and SME, besides the other SuperMAG in¬ 
dices, are primarily of the same nature, but with the Su¬ 
perMAG indices being enhanced by the higher number of 
ground based stations used to build the SuperMAG indices. 

More details about the SuperMAG initiative can be found in 
Gjerloev [2009]; Newell and Gjerloev [2011], and an expla¬ 
nation about data techniques and assimilation is reported 
by Gjerloev [2012]. Finally, the data are available from 
the SuperMAG websites http://supermag.jhuapl.edu/ and 
http://supermag.uib.no/. 

2.2. Determination of shock parameters and event 
analyses 

IP shocks during the period investigated here have 
been cataloged by several sources, such as the Havard- 
Smithsonian Center for Astrophysics (CfA) IP shock 
list compiled by Dr. J. C. Kasper located at 
http://www.cfa.harvard.edu/shocks/wi_data/ for WIND 
data, and http://www.cfa.harvard.edu/shocks/ac_master_data/ 
for ACE data. We also used a shock list com¬ 
piled by the ACE team available at http://www- 
ssg.sr.unh.edu/mag/ace/ACElists/obs_ list.html^shocks. An¬ 
other source used was the shock list with only ACE data 
from February 1998 to August 2008 published by Wang 
et al. [2010]. All these lists were merged to compile the 
shock list used here. We also used an automated search 
program to detect IP shock candidates in the raw data. 
After the shock was visually inspected, and if it satisfied 
the Rankine-Hugoniot conditions, the event was included in 
our list. Other sources were also consulted for comparison 
among several events in terms of solar wind conditions and 
IP shock parameters, such as calculated IP shock normal an¬ 
gles and speeds, when available [Berdichevsky et al., 2000; 
Russell et al., 2000; Zhou and Tsurutani, 2001; Pfech et al., 

2008; Wang et al., 2009; Richardson and Cane, 2010; Koval 
and Szabo, 2010; Zhang et al, 2012; Grygorov et al, 2014]. 

Once a shock was identified, solar wind data from WIND 
and ACE were inspected to provide the basis for IP shock 
parameter calculations. It is well known that IP shock nor¬ 
mal calculations are very sensitive to upstream and down¬ 
stream plasma parameters. Then, the highest quality avail¬ 
able spacecraft data were chosen for shock parameter deter¬ 
minations as described below. From a total of 461 identified 
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fast forward IP shocks, 272 were observed by ACE (59%), 
and 189 were observed by WIND (41%). 



Figure 1. The IP shock event seen by ACE on 23 
June 2000 at 1227 UT is an example of the methodol¬ 
ogy used for shock normal calculation and geomagnetic 
activity analysis. Plots (a)-(f) show jumps in B z and 
total magnetic held, in nT; thermal plasma pressure, in 
pPa; particle number density, in cm -3 ; shock speed, in 
km/s; and dynamic pressure DP (pv 2 ), in nPa. Plots 
(g) and (h) show SuperMAG data for the symmetric ring 
current SMR (similar to SYM-H), SME, and SML, all 
in nT. The maximum geomagnetic activity was recorded 
for both SME and SML approximately 2 hours after the 
shock impact. The time interval used to identify peaks 
in geomagnetic activity for all IP shocks was from ^30 
minutes to 2 hours after shock impacts. 


There have been a variety of shock normal determina¬ 
tion methods suggested since late 1960s. Some of the 
most commonly used methods using single spacecraft data 
are the magnetic and velocity coplanarity [Colburn and 
Sonett, 1966], and the mixed IMF and plasma data meth¬ 
ods [Abraham-Shrauner , 1972; Abraham-Shrauner and Yun, 
1976; Vihas and Scudder, 1986]. Although situations where 
data from more than one spacecraft are available give more 
reliable results [Burlaga et al., 1980; Russell et al., 1983a, b; 
Russell and Alexander , 1984; Thomsen, 1988; Russell et al, 
2000; Szabo, 2005; Koval and Szabo, 2010], we use the meth¬ 
ods based only on one spacecraft. Multiple spacecraft data 
usage would create an inconsistent data set in a large sta¬ 


tistical study, because the availability of more than one 
spacecraft data for shock normal determination is rare. The 
IP shock normals are then calculated using the methods of 
magnetic coplanarity (MC), velocity coplanarity (VC), and 
the three mixed data methods (MX1,MX2,MX3) found in 
Schwartz [1998]. Then, the average of the at least three 
closest results is calculated and registered as the chosen IP 
shock normal for each event. 

An example of an event analysis is shown in Figure 1. 
This shock event was seen by ACE on 23 June 2000. At 
1227 UT and (234, 36.6,-0.7 )Re GSE upstream of the Earth, 
ACE observed sharp jumps in magnetic field B z compo¬ 
nent, total magnetic field, plasma thermal pressure, parti¬ 
cle number density, plasma velocity, and dynamic pressure 
pv 2 (Figure 1 (a)-(f))- Approximately 55 minutes later, the 
shock impacted the magnetopause, the magnetosphere was 
compressed by the shock and an SSC was detected by Super¬ 
MAG geomagnetic stations, as can be seen in Figure 1(g) for 
the SuperMAG symmetric ring current index SMR [Newell 
and Gjerloev, 2012], the SuperMAG index similar to the well 
known SYM-H index. Increases in the SuperMAG indices 
SME and SML followed the IP shock approximately 2 hours 
after shock impact, reaching a maximum of about 1500 nT 
for SME and a minimum of about -1000 nT for SML. The 
calculated shock normal of this event is (-0.785,0.153,-0.600), 
with 9 Xn ^140°, shock speed of 553.2 km/s, and fast magne- 
tosonic Mach number 2.60. Using these results, and assum¬ 
ing the estimated position of the magnetopause previous to 
the shock impact as 10 Re as suggested by Zhou and Tsu- 
rutani [1999], the calculated time travel is ~55 minutes, in 
agreement with observations, which validates our method. 
To complete the shock property analysis, the compression 
ratio (the ratio of downstream to upstream plasma density) 
was 2.62, and the fast magnetosonic Mach number was 2.60. 

3. Statistical results 

3.1. Solar wind and shock parameters 

Figure 2 shows the yearly IP shock number (gray bars) 
and the monthly sunspot number (SSN, solid lines) plot¬ 
ted in the time range from 1995 to 2013. This time period 
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Figure 2. Yearly IP shock number (gray bars) plotted 
against the SIDC monthly sunspot number (solid line). 
WIND and ACE data were used to identify all IP shock 
events. A strong correlation can be seen. The maximum 
yearly IP shock number occurred in the year 2000 (50 
events), in the solar maximum of the solar cycle 23. Due 
to the unusually low sunspot number in the maximum of 
the current solar cycle, only 25 events were observed in 
2013, and not many more are expected to be identified 
in the 2014 and 2015 WIND and ACE data [Smith et al., 
2014], 
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includes the whole solar cycle 23, which ranged from May 
1996 to January 2008. A correlation between both numbers 
is clear. During the ascending phase of the solar cycle 23, the 
number of IP shocks increased with the SSN. Correlations 
between the number of SSCs and the SSN in different solar 
cycles have been reported by earlier works [Chao and Lep- 
ping, 1974; Hundhaunsen, 1979; Smith, 1983; Smith et al., 
1986; Rastogi, 1999]. Since 75% of SSCs are caused by IP 
shocks [Wang et al, 2006], these arguments are considered 
to be very similar. Such correlations may be a result of the 
increasing number of shocks driven by ICMEs [Jian et al., 
2006a], since the number of CIR-driven shocks did not show 
any correlation with the solar cycle phase [Jian et al., 2006b], 
Then, during the declining phase of the solar cycle 23, the 
number of IP shocks decreases with the SSN, since the num¬ 
ber of ICMEs related to shocks decreases. The number of 
fast forward shocks is higher in the maximum of the solar 
cycle 23 (year 2000, 50 events) than in the solar minimum 
of the same solar cycle (years 1995-1996, 31 events) [Echer 
et al, 2003]. Due to the unusual low SSN of the current solar 
cycle maximum, barely more than 25 shocks are expected 
to be found in the WIND and ACE data for 2014 and even 
2015 [Smith et al., 2014], 

A statistical analysis of solar wind and IP shock param¬ 
eters is shown in Figure 3(a-f). Figure 3(a) shows 9 X „, the 
angle between the shock normal vector and the Sun-Earth 
line. Angles close to 180° indicate IP shocks whose shock 



Figure 3. Statistical results of the full list with 461 
IP shocks. Figure 3(a) shows the distribution of the im¬ 
pact angle between the shock normal and the Sun-Earth 
line. Angles close to 180° represent almost frontal shocks. 
The clock angle ip yrl distribution on the GSE YZ plane is 
shown in Figure 3(b). Angles in the ranges 0° < (p y „ < 
45°, 135° < ifi Vn < 225°, and 315° < ip Vn < 360° indi¬ 
cate that the shock normals were close to the equatorial 
plane. Figure 3(c) shows the distribution of 9b„, the an¬ 
gle between the upstream magnetic held vector and the 
shock normal. Represented in Figure 3(d), is the shock 
speed (in km/s) distribution, in relation to the spacecraft 
frame of reference. Figure 3(e) shows the distribution of 
the compression ratio, the ratio of the downstream to up¬ 
stream plasma densities. Finally, the fast magnetosonic 
Mach number distribution is shown by Figure 3(f). 


normals were almost parallel to the Sun-Earth line, or al¬ 
most frontal shocks. IP shocks with angles close to 90° rep¬ 
resent inclined shocks. In our list, 363 (78.57%) cases had 
shocks with 9 Xn > 135°. The distribution of the clock an¬ 
gle ip Vn is shown in Figure 3(b). Shock normals with 0° 
< <Py n < 45°, 135° < ip y „ < 225°, and 315° < ip Vn < 
360° indicate that the shock normal was close to the equa¬ 
torial plane. These conditions were satisfied by 276 events, 
or 59.74%. Figure 3(c) shows the obliquity 9s n , the angle 
between the shock normal and the upstream magnetic held 
vector. In our data set, 354 cases showed larger than 
45°, and most of the shocks in this category might have 
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Figure 4. SML jumps, in nT, triggered by IP shock im¬ 
pacts are plotted as a function of the shock speed v s , in 
km/s. The events were binned in three different groups 
in terms of the shock orientation in relation to the Sun- 
Earth line: Figure 4(a), 120° < 9 Xn < 140° (highly in¬ 
clined shocks); Figure 4(b), 140° < 9 Xn < 160° (inclined 
shocks); and Figure 4(c) 160° < 9 Xn < 180° (almost 
frontal shocks). The IP shocks are more geoeffective for 
strong (high speed) and almost frontal shocks (large 9 Xn ). 
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been driven by ICMEs [Richardson and Cane, 2010]. The 
shock speed distribution is shown in Figure 3(d). The av¬ 
erage shock speed is 467 km/s, and it tends to be higher in 
solar maximum, and lower in solar minimum, as already re¬ 
ported by Berdichevsky et al. [2000] and Echer et al. [2003] 
with data partially in the same time period. The percent¬ 
age of shocks above the average speed is 40.13%, or 185 
events. The compression ratio, the ratio of the downstream 
to upstream plasma densities, can be seen in Figure 3(e). 
As reported before [Berdichevsky et al, 2000], most shocks 
have their compression ratios between 1.2 and 2.0, which 
happened to 251 of our cases (54.44%). Our compression 
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ratio average is 2.07. Although the theoretical limit for the 
compression ratio is 4 [Richter et al, 1985], which is de¬ 
rived for perpendicular shocks, this value was exceeded in 
11 cases (2.38%), and most of them took place slightly be¬ 
fore and after the solar maximum (year 2000). Echer et al. 
[2003] argued that such cases can happen for some shocks 
in a data set in which shock obliquities range from almost 
parallel to almost perpendicular shocks. Finally, the fast 
magnetosonic Mach number distribution is shown in Figure 
3(f). The average of Ms is 2.15, and it is clear that most 
shocks have Ms between 1.0 and 2.0 [Kallenrode, 2003]. The 
number of shocks with Ms above the average is 166 (36.00%). 
However, some shocks have Ms less than one, which can be 
an indication that such events were not shocks because the 
shock waves could not steepen, even though they could show 
some shock-like behavior [Kennel et al, 1985]. These events 
were not included in our statistical analysis. Therefore, as a 
consequence of this analysis, it is possible to conclude that 
the interplanetary space is dominated by weak IP shocks. 
The agreement of our results with other works validates our 
statistical analysis, in particular the shock normal determi¬ 
nation methods used in this work. 

3.2. Geomagnetic activity 

In this section we investigate the geoeffectiveness of IP 
shocks by correlating the shock parameters with the Su- 
perMAG SML index as a geomagnetic activity indicator. 
Changes in this index, ASML, in nT, are recorded for each 
event from ~30 minutes to two hours after shock impact. If 
the IP shock is followed by any other solar wind structure, 
only the first peak in the data is considered. We chose this 
time frame because some inclined shocks take a long time 
to sweep over the magnetosphere when they are inclined in 
relation to the Sun-Earth line [Takeuchi et al., 2002; Guo 
et al., 2005; Wang et al., 2006; Oliveira and Raeder, 2014]. 
We used SuperMAG data up to 2013 because the 2014 Su- 
perMAG data were not yet available. 

Figure 4 shows jumps in SML, in nT, measured by Su¬ 
perMAG ground stations plotted against the shock speed, 
in km/s. Since we consider two parameters, shock speed 
and impact angle, all the data were binned in three dif¬ 
ferent groups in terms of the shock normal impact angle 
8 Xn . Here, the impact angle is held and the shock speed 
varies. Figure 4(a) shows highly inclined shocks: 120° 

< 9x n < 140°; Figure 4(b) represents moderately inclined 
shocks, 140° < 8 Xn < 160°; and almost frontal shocks, 160° 

< 9 Xn < 180°, can be found in Figure 4(c). In Figure 4(a), 
most shocks produce little geomagnetic activity (ASML < 
500 nT), and in such cases most shocks had v s < 450 km/s. 
This is expected for weak and highly inclined shocks. For 
some stronger, but highly inclined shocks, the resulting ac¬ 
tivity is slightly larger, but just a few such shocks in this case 
were identified in the data. The linear regression analysis 
gives a correlation coefficient of R = 0.34. In the intermedi¬ 
ate case, i.e., the case of shocks with moderate inclination, 
most shocks produced ASML > 500 nT. In this case, there 
is a stronger correlation. We attribute the correlation coef¬ 
ficient of R = 0.41 to the fact that most shocks with v s < 
450 km/s triggered small jumps in SML (ASML < 500 nT). 


Figure 5. SML jumps, in nT, triggered by IP shock 
impacts are plotted as a function of the shock impact 
angle 0 Xn , in km/s. The events were binned in three dif¬ 
ferent groups in terms of the shock speed: Figure 5(a), 
350 < v s < 450 km/s (weak shocks); Figure 5(b), 450 
< v s < 550 km/s (moderate shocks); and Figure 5(c) 
v 3 > 550 km/s (strong shocks). The shocks are more geo¬ 
effective for almost frontal (large 9 Xn ) and strong (high 
speed) IP shocks. 


Table 1. Summary of the results obtained for the shock 
speed, shock impact angle, and ASML correlation analyses. 


Fixed 

impact angle 6 X 

n , changed shock speed v s 

category 

R 

highly inclined moderately inclined 
0.42 0.47 

almost frontal 
0.67 

Fixed shock speed v s , 

changed impact angle 6 Xri 

category 

R 

weak 

0.37 

moderate 

0.48 

strong 

0.78 
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For the cases in which v s > 450 km/s, ASML showed better 
correlations, but just a few with ASML > 1000 nT. In the 
more extreme case, namely the case in which the IP shocks 
were almost frontal, the correlation coefficient is R = 0.60. 
In this case, approximately half of the shocks with v s < 
450 km/s did not show large jumps in SML. Most shocks 
triggered ASML > 500 nT, and almost all cases in which 
ASML > 1000 nT had v s larger than 450 km/s. Therefore, 
by inspecting all plots, it is clear that the IP shock geo¬ 
effectiveness increases with both shock strength and shock 
impact angle. Table 1 summarizes the results obtained in 
all categories in this case. 

The opposite analysis is shown in Figure 5, i.e., the shock 
speed is held and the impact angle varies. There, ASML 
is plotted against 9 Xn , an d the data are binned in three 
different categories related to the shock strength (or shock 
speed). Figure 5(a) shows the weak shocks, 300 < v a < 450 
km/s; 5(b) moderate shocks, 450 < v a < 550 km/s; and 
5(c) strong shocks, v s > 550 km/s. Figure 5(a) shows the 
largest number of small ASML (ASML < 500 nT), even for 
shocks with shock normals almost parallel to the Sun-Earth 
line. The correlation coefficient in this case is R = 0.36. A 
clearer ASML-6G n correlation is evident in the intermediate 
case, where R = 0.51, and most shock events have ASML 
> 500 nT and 9 Xn > 135°. All shocks with ASML > 1000 
nT had impact angles larger than 140°. In the category 
of strong shocks, only a few shocks triggered geomagnetic 
activity with ASML < 500 nT, most of them being highly 
inclined shocks in which 9 Xn < 150°. Shocks with high 
geoeffectiveness, or ASML > 1000 nT, were almost frontal 
shocks with 9 Xn > 150° (only one event had 9 Xn slightly less 
than 150° in this case). The highest correlation coefficient, 
R = 0.70, occurs for IP shocks in this category. Table 1 sum¬ 
marizes the results obtained in all cases in this correlation 
analysis. 

Thus, strong shocks are generally much more geoeffective 
than weak shocks, and the geoeffectiveness increases if the 
IP shock impacts more frontally the Earth’s magnetosphere. 
These results have already been shown by Wang et al. [2006] 
for the SSC rise-time and Oliveira and Raeder [2013, 2014] 
in global MHD simulations. 

4. Summary and Conclusions 

We investigated WIND and ACE solar wind data at 1 AU 
to compile a list of fast forward interplanetary (IP) shocks. 
We studied the geoeffectiveness triggered by the IP shock 
impacts, as measured by the jumps in the SuperMAG SML 
index, and how it relates to the shock speed (strength) and 
the shock inclination angles. Our main results are summa¬ 
rized below: 

1. We provide the community with a fast forward IP 
shock list with events from January 1995 to December 2013, 
covering the whole solar cycle 23 and half of the current 
solar cycle. 

2. The number of yearly IP shocks correlates closely with 
the monthly sunspot number. The highest number of fast 
forward IP shocks was found in the year 2000, the solar max¬ 
imum of the solar cycle 23. As expected, the number of IP 
shocks is smaller in the maximum of the current solar cycle 
due to the unusual low number of sunspots occurring in this 
period. 

3. The majority of the events (76%) are almost per¬ 
pendicular shocks, with 9s n > 45°. Most shocks (78%) 
have their shock normals close to the Sun-Earth line, or 
9 Xn > 135°. Also, less than half of the shocks (40%) have 
their speeds above the average of about 450 km/s, and 
shocks with the supermagnetosonic Mach number greater 
than the average 2.1 was 36%. These results indicate that 
the heliosphere at 1 AU is dominated by weak interplanetary 
shocks. 


4. Strong (high speed) shocks are more geoeffective than 
weak shocks (low speed). The correlation is clearer when 
shocks are grouped in categories related to their strength 
and then investigated in terms of their shock impact angles. 
The largest correlation occurs (R = 0.70) when we fixed the 
IP shock speed and changed the IP shock impact angles. 
Thus, the IP shock impact angle is just as important deter¬ 
mining their geoeffectiveness as their strength. This result 
was predicted by Oliveira and Raeder [2014]. 
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